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Abstract 

In this paper, we consider the output synchronization problem for a 
network of heterogeneous diffusively-coupled nonlinear agents. Specif- 
ically, we show how the (non-identical) agents can be controlled in 
such a way that their outputs asymptotically track the output of a 
prescribed nonlinear exosystem. The problem is solved in two steps. 
In the first step, the problem of achieving consensus among (identical) 
nonlinear reference generators is addressed. In this respect, it is shown 
how the techniques recently developed to solve the consensus problem 
among linear agents can be extended to agents modeled by nonlinear 
d-dimcnsional differential equations, under the assumption that the 
communication graph is connected. In the second step, the theory of 
nonlinear output regulation is applied in a decentralized control mode, 
to force the output of each agent of the network to robustly track the 
(synchronized) output of a local reference model. 



1 Introduction 

The problem of achieving consensus (among states or outputs) in a (homoge- 
neous or heterogenous) network of systems has attracted a major attention 
in the last decade. An exhaustive coverage of the literature, which is beyond 
the scope of the present paper, can be found, e.g. in the recent dissertation 
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|16j and in all references cited therein. We limit ourselves to mention that 
the case of a network of linear systems connected through a time-invariant 
graph has been fully addressed in the papers [13], [Hj, [T7], [6J, while the 
analysis of a network of linear systems connected through a time-varying 
graph reposes on a fundamental convergence result established in [ID], [11] . 
Major results concerning the consensus problem in a network of nonlinear 
systems can be found in [5] , p2] , [7] , [I] , [IS] . The purpose of this paper is 
to present further contributions to the problem of output synchronization 
in a network of heterogeneous diffusively-coupled nonlinear agents. 

As shown in [17] for linear systems and in [16] for nonlinear systems, if 
the outputs of the agents of a heterogenous network achieve consensus on 
a nontrivial trajectory, the trajectory in question is necessarily the output 
of some autonomous (linear or nonlinear, depending on the case) system. 
This is the equivalent, in the context of the consensus problem, of the cel- 
ebrated internal model principle of control theory |3]. Motivated by this, 
we consider in what follows the problem of controlling a set of networked 
(non-identical) nonlinear agents in such a way that their outputs asymptot- 
ically track the output of a prescribed nonlinear exosystem. The problem is 
solved in two steps. In the first step, a network of N identical copies of the 
given nonlinear exosystem is considered, the k-th of which is to be seen as 
"local reference generator" for the k-th agent, and it is shown how certain 
"coupling gains" can be chosen in such a way that these local generators 
synchronize on a common consensus trajectory. To this end, we extend ex- 
isting techniques (see [13] . [14] ) recently proposed for the synchronization of 
a homogeneous network of linear systems exchanging information through 
a connected (time-invariant) communication graph. The arguments used in 
this part are inspired by the literature on high-gain stabilization of nonlin- 
ear systems by output feedback, and specifically by the design of high-gain 
observers (see e.g. [1]). In the second step, the problem of controlling the 
individual agent in such a way that its output tracks a reference output 
generated by the "local exosystem" is addressed as a "classical" problem of 
nonlinear output regulation. In this respect, it is shown how the theory of 
nonlinear output regulation proposed in [9] (see also [2] , [8] ) , can be success- 
fully adopted to design robust internal model-based local regulators if the 
dynamics of the local agents fulfill a (weak) minimum-phase assumption. 
The result presented in the paper can be seen as a kind of "separation prin- 
ciple" , in which the tools used to design local regulators having the internal 
model property with respect to a local nonlinear exosystem in steady state, 
and the tools for synchronizing a set of networked nonlinear homogenous 
exosystems to reach a common steady state, can be combined to achieve 



consensus of the outputs of heterogenous networked nonlinear systems. 

The paper is organized as follow. In the next Section the problem is pre- 
cisely formulated and the structure of the controllers is specified. Section III 
presents the solution for the first of the steps mentioned above, namely con- 
sensus in a network of diffusively-coupled identical nonlinear systems. The 
problem of reaching a consensus between heterogenous systems by means of 
nonlinear internal model-based regulator is addressed in Section IV, while 
Section V presents some simulation results, concerning the theory presented 
in Section III. 

2 Problem statement 
2.1 Communication graphs. 

In what follows, the communication between individual systems (agents) is 
encoded by a time-invariant communication graph. The latter is a triplet 
G = {V, E, A} in which: 

• V is a set of N vertices V = {v±,V2, ■ ■ ■ ,vn}, one for each of the ./V 
agents in the set. 

• £ C V x V is a set of egdes that models the interconnection between 
nodes, according to the following convention: (vk,Vj) belongs to £ if 
there is a flow of information from node j to node k. 

• the flow of information from node j to node k is weighted by the (k,j)- 
th entry a^j of the adiacency matrix A G R NxN . 

It is assumed that there are no self-loops, i.e. that (vk,Vk) £ £. The set 
of neighbors of node Vk is the set Nk = {vj G V : a^j / 0}. A path from node 
Vj to node Vk is a sequence of r distinct nodes {v^ , . . . , vg r } with vg 1 = Vj 
and V£ r = Vk such that (vi + ±,Vi) G £. A graph Q is said to be connected if 
there is a node v such that, for any other node Vk G V \ {v}, there is a path 
from v to Vk- 

In what follows, we will consider cases in which the information available 
for control purpose at the k-th agent at time t has the form 

N 

«* = Z)ajy(<i(*)-Cfc(*)) k = l,...,N (1) 



in which Q, for i = 1, . . . , N, is a measurement taken at agent i. Letting L 
denote the so-called matrix Laplacian matrix of the graph, defined by 




-a kj (t) for k^j 



N 



^4,0(t) k = i,...,N, 




By definition, the diagonal entries of L are non-negative, the off-diagonal 
elements are non-positive and, for each row, the sum of all elements on this 
row is zero. As a consequence, the all-ones N- vector ljv = col(l, 1, . . . , 1) 
is an eigenvector of L, associated with the eigenvalue A = . Let the other 
(possibly nontrivial) N — 1 eigenvalues of L be denoted as X^(L), . . . , Xn(L). 

Theorem 1 A time-invariant graph is connected if and only if its Laplacian 
matrix L has only one trivial eigenvalue Ai = and all other eigenvalues 
X2(L), . . . , Atv(L) have positive real parts. 

2.2 Problem formulation 

We consider in what follows the problem of inducing consensus between the 
outputs of N non-identical nonlinear systems, which exchange information 
through a communication graph Q. The control system is decentralized, 
i.e. there is no leader sending information to each individual system, but 
rather each system exchanges information only with a set of neighboring 
systems, the information in question concerning only the relative values of 
the respective controlled outputs. The N nonlinear agents are described 



k = 1, . . . , N, where u k and y k are the local control input and output, with 
the inputs u& that must be designed in such a way that the outputs of 
the iV systems asymptotically reach consensus on a nontrivial common tra- 
jectory y*(t). Each agent is controlled by a local output-feedback controller 
of the form 



by 





(3) 



Uk = ik{Ck,Vk,v k ) u k eR 
Cfc = Pk{£k,yk) Ck g R p 



(4) 



in which ^ and i>k are outputs and inputs that characterize the exchange 
of relative information between individual (controlled) agents, which takes 
the form (pQ). 

In general terms, the control problem can be formulated as follows. Let 
Xk G M rafc , k = 1,...,N, be fixed compact set of admissible conditions for 
([3]). The problem is to find N local controllers of the form @, exchanging 
information as in ([1]), and compact sets £ R nfc , k = 1, . . . ,N, of admis- 
sible initial conditions for all such controllers, so that the positive orbit of 
the set of all admissible initial conditions is bounded and output consensus 
is reached, i.e. for each admissible initial condition (xj.(0), £fc(0)) € x 
k = 1,...,N, there is a function y* : R — > R d such that 

lim \y k (t) - y*(t)\ = for all k = 1, . . . , N , 

t— >oo 

uniformly in the initial conditions. 

With the results of [16] in mind, we expect that the consensus trajectory 
y*(t) can be thought of as generated by a nonlinear autonomous system, 
which could be modeled as an ordinary differential equation of order d 

2 /*( d ) = 0(y*,y*( 1 ),...,y*( d - 1 )) (5) 

or in the equivalent state-space form of a d-dimensional system with output 

w = s(w) , w G R d , x 

y* = ${w) [ ' 

in which 

s{w) = Sw + B<t>(w) , #{w) = Cw (7) 

and (S,B,C) is a triplet of matrices in prime form. Since we are seeking 
nontrivial consensus trajectories, in what follows we will consider the case in 
which © possesses a nontrivial compact invariant set W. Moreover, we will 
assume that the function </>(•) is globally Lipschitz. In presence of systems of 
the form (|7|) in which the </>(■) is only locally Lipschitz, this assumption can 
be always enforced by properly modifying the function outside the compact 
set W by using appropriate extension theorems. 

2.3 Structure of local controllers and communication proto- 
col 

Bearing in mind the possibility of modeling all solutions of ([5]) as outputs of 
the autonomous system ([S])-©, in what follows, we consider for the local 



controllers (j3|) a structure of the form 



N 

w k = s(w k ) + Ky^ j a kj ('&(wj) -ti(w k )) 

(8) 

Vk = <Pk(Vk,ek) 
Uk = lk{Vk,e k ) 

in which 

ek = Vk~ d(w k ) • (9) 

It is readily seen that this structure consists of a set of N local reference 
generators 

w k = s(w k ) + Kv k , , 

y k ^ = d{w k ) , [W) 

coupled via 

N 

each one of which provides a reference y k ei to be tracked by a local regulator 

Vk = fk(Vk,ek) 
Uk = -Yk(Vk,ek) 

driven by the local tracking error 

e-k = Vk-Vk ■ 

This control structure enables us to solve the problem in two stages. In 
the first stage, the design parameter K is chosen in such a way as to induce 
consensus among the local generators fjlOj) . In the second stage, the local 
regulators are designed in such a way that each of the outputs yk tracks 
its own reference y k ret . It goes without saying that in the second step will 
ought to be able to use - off the shelf - a large amount of existing results 
about the design of output regulators for nonlinear systems in the presence 
of exogenous signals generated by a nonlinear exosystem. 

In this framework, we address first the problem of achieving consensus 
among the N local generators ()10p . Similar problems have received a large 
attention in recent literature (an exhausting covering of all such literature 
is beyond the scope of this paper, excellent surveys can be found in [16], 
[13] . [2]) but, to the best of our knowledge, a solution to the problem in 
the general terms considered here, i.e. for a network of N nonlinear systems 



of dimension d > 1, with information exchange in terms of relative values 
of 1-dimensional outputs (other than relative values of their d-dimensional 
states), has not been proposed yet. We will address this problem in the 
following section. 

3 Achieving consensus in a homogeneous network 
of nonlinear systems 

As anticipated, we consider in what follows the problem of achieving state 
consensus in a network of N identical nonlinear systems of the form ()10p 
coupled as in (fTT|) . in which s(w) and $(w) are the map and the function 
defined in ((7j). 

The problem will be solved under the following assumptions. 
Assumption 1 The graph Q is connected. 

Assumption 2 There exists a compact set W cM. d invariant for (0|) such 
that the system 

w = Sw + Bcj}(w) + v 

is input-to-state stable with respect to v relative to W , namely there exist a 
class-JCC function /?(•,•) and a class-JC function 7(-) such that 

\\w(t,w)\\ w < max{P(\\w\\ w ,t), j( sup ||«(r)||)}. 

re[0,t) 

To the purpose of inducing consensus in the network ([10]) — fpTT]) . we choose 
the vector K in (fTUl) as 

K = D g K , 

where D g = diag ( g g 2 ... g d ) , with g a design parameter and Kq a 
vector to be designed. 

Due to Assumption 1, it is known that the Laplacian matrix L has only 
one trivial eigenvalue and all remaining eigenvalues have positive real part. 
Hence there exists a fx > such that 

Re[Aj(L)]>/i i = 2,...,N. (12) 

With this in mind, let T G W NxN be defined as 

T = f i o r _^ 

V J-JV-1 1 N-l J 



and note that 

L = T- 1 LT=( n ° L T 21 

V UAT-ixl -^22 

in which the eigenvalues of L22 coincide with X^L), . . . , X^(L). Then, the 
following result holds. 

Lemma 1 Let P be the unique positive definite symmetric solution of the 
algebraic Riccati equation 

SP + PS T - 2fiPC T CP + al = 

with a > 0. Take Kq as 

K = PC T . (13) 

Then, the matrix 

[(In-! ® 5) - (L 22 ® K C)] 

is Hurwitz. < 

The proof of this Lemma can be found in |14j or in [16]. Using this, we 
can now proceed with the proof of the main result of this Section. 

Proposition 1 Suppose Assumptions^ and\^ hold. Consider the network 
of N coupled systems 

N 

w k = Sw k + B(j)(w k ) + D g K ^2 akj(Cuij - Cw k ) (14) 

i=i 



with k = 1, . . . , N. Let Kq be chosen as in H3\). Then, there exists a number 
g* > such that, for all g > g* , the compact invariant set 

W = {(w!,W2, . ■ ■ ,wn) eWxWx---xW: . , 

W! = W2 = ■ ■ ■ = VJNJ 

is globally asymptotically stable. < 



Corollary 1 Let the hypotheses of the previous Proposition hold and let Kq 
be chosen as in \13\) . There is a number g* > such that, if g > g* , the 
states of the N systems JT^P reach consensus, i.e. for every w k (0) £ 
k = 1, . . . , N, there is a function w* : R — > M. d such that 



lim \w k (t) -w*(t)\ = for all k = 1,... ,N . < 

t— ¥00 



Proof. By the definition of Laplacian, the fc-th controlled agent of the 
network (1141) can be written as 



N 

w k = Sw k + Bcj)(w k ) - DgKo^lkjCwj . 

3=1 

Thus, setting w = col(u/i, . . . , wn) entire set of iV controlled agents can be 
rewritten as 

W = [(Ijv ® S) - (L <g> D g K C)} w+(I N ® 5)$(w) 

where 

$( w ) = col(<f>(wx), 4>(w N )) . 

Consider the change of variables w = (T _1 ® /d)w, in which T is the 
matrix introduced above. The system in the new coordinates read as 

ib = (T- 1 ®I d )[(I N ®S)-(L®D g K C)](T®I d )w 
+ (T" 1 ® I d )(I N ® B)$((T ® J d )w) 



(IiV ® 5) - (L ® Z) 9 iTo<?) 



w 



+ (T _1 ® B)$((T ® 7(i)w) . 



Observing that 



w = col(wi, — u>i, . . . , wat — u>i) 
set z k = Wfc — wi, for k = 2, 3, . . . , N, and 

z = col(z 2 ,^3, ■ ■ ■ ,2jv) , 

yielding w = col(wi,z). Then, it is readily seen that the system above 
exhibits a triangular structure of the form 

wi = Swi + B(f>(wi) - (L 12 ® D g K Q C)z 
z = [(I N „ 1 ®S)-(L 22 ®D g K C)]z + A<S>(wi,z,) 



where 



A$(u7i,z) = (Jjv-i ® -B) 



/ + z 2 ) - 0(^i) \ 

+ z 3 ) - 4>{wi) 

\ (j)(w 1 + z N ) - (p(wi) J 



Note that A<&(w\, z) is globally Lipschitz in z uniformly in w\ and A<£(u>i, 0) = 
for all 6 M n . Consider now the rescaled state variable 

C = (/at-i®1VV 
By using the definition of S, C and B, it follows that 

wi = Sw x + B(p( Wl ) - (L12 ® D g K C)(I N -i <g> D 9 )C 
C = 5 [(/jv-i 5) - (L 22 K C)] C 

+ -,A$(u>i,(J^_i®Z> ff )C). 

It is known from the Lemma above that the proposed choice of Kq guarantees 
that the matrix [(Ijv-i <8> S) — (L 22 <8> KqC)] is Hurwitz. As a consequence 
standard high-gain arguments lead to the conclusion that if g is chosen 
sufficiently large, the equilibrium Q = of the lower subsystem is asymp- 
totically stable, uniformly in wi, and actually with a quadratic Lyapunov 
function that is independent of w\. The ISS property in Assumption 2 thus 
guarantees that w\ converges to the invariant set W. Since Wk = w\ + Zk, 
k = 2, . . . , N, and zt — > as t — > oo, the result follows. < 



4 Achieving consensus in a heterogenous network 
of nonlinear systems 

We proceed now with the second step of the design, i.e. we design local 
regulators for each agent. In what follows, we assume that the vector Kq and 
the values of g have been fixed and such that the conclusion of Proposition 
[T] holds, i.e. such that the set (fT5j) is globally asymptotically stable for (fT4|) . 

In what follows, we assume that each individual agent has a well defined 
relative degree r between input and output and possess a globally 
defined normal form. To streamline the exposition, we consider the special 
case in which r = 1. The case of higher relative degree only entails heavier 
notational complexity and no conceptual differences. Thus we assume that 
the individual agent is modeled by equations of the form 

Zk = fk(Zk,Vk) Qgs 

Vk = qk{zk,Uk) + bk(zk,Vk)uk , 

where z/% € M"* -1 and where bk(zk,yk), which is the high-frequency gain of 
the k-th agent, is bounded away from zero. In particular, we assume that 



there exists a b k > such that b k (z k ,y k ) > b k for all (z k ,y k ) G M nfc ^ 1 x R 
and for all k = 1, . . . , N. As anticipated, with this system we associate a 
local tracking error of the form 



The problem is to design a local regulator to the purpose of steering e k to 
zero. 

In this respect, it should be borne in mind that Wk is a "portion" of 
the state of the coupled system (|14|) and hence the entire dynamics of the 
latter should be taken into account in the analysis. To this end, invoking the 
arguments used in the proof of Proposition [H observe that, for any k, there 
exists a linear change coordinates in (|14p . by means of which this system 
(the entire set of N networked local reference generators) can be changed 
into a system modeled by equations of the form 



in which ip(0,w k ) = 0. By assumption, the lower sub-subsystem is input- 
to-state stable, with respect to the input (k, to the set W. Moreover, as 
observed in the proof of Proposition [H the equilibrium ( k = of the upper 
subsystem is globally exponentially stable. As a consequence, the set 



is a globally asymptotically stable compact invariant set of (fT7|) . 

In view of this, we can represent the aggregate of (fTHj) and of (fTTl) as a 
standard "exosystem-plant" interconnection 



ek = Vk- Cw k . 



w k = s(w k ) + T fc Cfc 



(17) 



{(0 



■k,w k ) 6 R^" 1 ^ xR d :( k = 0,w k £ W} 



Ck = 
w k = 



Vk = 



e-k = 



Zk = 



ip(Ck,w k ) 
s(w k ) + T fc Cfc 
fk(z k ,yk) 

Qk(zk,yk) + b k (z k ,y k )u k 
y k - Cw k ■ 



(18) 



As usual, we change variables replacing y k by e k and obtain 



Ck 

w k 



Tp{(k,w k ) 
s(w k ) + T fc Cfc 
fk(z k , Cw k + e k ) 

q k (z k , Cw k + e k ) - C[s(w k ) + T k ( k ] 



(19) 



+ b k {z k ,Cw k + e k )u k . 



This system is ready for the design (under appropriate hypothesis) of a local 
regulator of the form 

Vk = fkiVk) + G k v k 

Uk = lk{rik)+Vk (20) 
Vk = K k (e k ) 

according to the procedures suggested in [2] or in [9] . The basic assumption 
needed to make the design possible is that the zero dynamics of (|19p namely, 
those of 

Ck = Tp{(k,w k ) 

Wk = s(w k ) + TkCk (21) 
Zk = fk(z k ,Cw k ) 

possess a compact invariant set which is asymptotically stable with a domain 
of attraction that contains the prescribed set of initial conditions. To make 
this assumption precise, let Wk be the set of admissible conditions of Wk, 
let Sk be the set of admissible initial conditions of Ck and Zk the set of 
admissible initial conditions of Zk- Then, the assumption is question can be 
stated as follows. 

Assumption 3 There exist a (possibly set-valued) map 7Tfc : w/~ G W t— > 
7Tk( w k) C M"* -1 such that the set 

■Ak = {{(k,Wk,z k ) ■■ Ck = 0, w k G W, z k G -K k (w k )} 

is an asymptotically stable invariant set for $21\) with a domain of attraction 
containing Sk x Wk x Zk- 

We note that this assumption is the natural formulation, in the current 
framework of a networked system, of the (weak) minimum-phase assumption 
that one would assume in solving a problem of output regulation for the fc-th 
agent if high-gain arguments were to be used for stabilization purposes. 

We proceed now with the design the functions (ipk(-),^fk(-),Gk) in (|20j) . 
whose key properties are captured in the following definition, taken from [8]. 

Definition (Asymptotic internal model property). The triplet (yfc(-), 7fc(-)j Gk) 
is said to have the asymptotic internal model property if there exists a C 1 
map T k : R d x W 1 ^ 1 -)■ M m fe such that the following holds: 

x We use the notation gr(7r fc ) := {(w k , z k ) : w k £ W, z k € Tv k (w k )} 



(i) for all (w k ,Zk) € gr(7T fc ) 



(ii) the set 



OWk oz k 

Cs(w k ) - q k (z k ,Cw k ) 
b k (w k ,Cw k ) 



7k(T(w k ,Z k )) 



Sk = {((k,w k ,z k ,r] k ) : 

C k = 0, (w k ,z k ) G gr(7r fc ), r] k = r(w k , z k )} 

is locally asymptotically stable for the system 



Ck 

Wk 

Zk 

Vk 



^(Ck,w k ) 
s(w k ) + T fc Cfc 
f k (z k ,Cw k ) 

<Pk(Vk) ~ G k 7fc(%) 



+ 



q k (z k , Cw k ) - C[s(w k ) + T fc Cfc 



b k (w k ,Cw k ) 



with a domain of attraction containing Sfc x W k x Z k x M k , where M k is the 
compact set of initial conditions of (|20p . < 

If a triplet with the asymptotic internal model property can be designed 
then the problem of steering the regulation error e k of the k-th agent to zero 
is solved as claimed by the following theorem proved in [9]. 



Theorem 2 Let S k C R( Ar ~ 1 ) (i ; W k C R d , Z k C R"* -1 , E k C R and 
Mfc C W rik be compact sets of initial conditions for the closed-loop system 
[TPj) . IW\) . Let the triplet {ip k {-), ~y k (-), G k ) be designed so that it has the 
asymptotic internal model property. Then there exists a continuous function 
K k : K — > M such that the trajectories of the closed-loop system originating 
form S k x W k x Z k x E k x M k are bounded and lim^oo e k {t) = uniformly 
in the initial conditions. 



A triplet having the internal model property can be always designed as 
detailed in the next result coming from a slight adaptation of the results 
presented in [9]. 



Proposition 2 Let m k > 2(d + n k — 1) + 2. Then there exists a < 
and, for almost all possible choice of controllable pairs (F k ,G k ) € W nkXnik x 
l™ tXl such that Rea(F k ) < \* k , there exists a continuous ^ k : R mfc — > R, 
such that the triplet {tp k (-), 7fc(-)> G k) with <p k (r] k ) = F k r] k + G k j k (r] k ) has 
the asymptotic internal model property. 

The previous results, although conceptually interesting, is not construc- 
tive in the design of the function j k . As shown in [2], it turns out that a 
constructive design procedure can be given if an extra assumption is invoked. 
In particular, assume that there exists a m k > and a locally Lipschitz func- 
tion g k : R mft — y R with the property that, for all (w k (0), z k (0)) € gr(ir k ), 
the solution w k (t), z k (t) of 

w k = s(w k ) 
z k = fk{zk,Cw k ) 

is such that the function 

^ _ Cs(w k (t)) - q k (z k (t),Cw k (t)) 
b k (w k (t),Cw k (t)) 

satisfies 

p {mk \t) = Q k {p{t), P {l) {t),...,p^- l \t)) Vt€R. 

If this assumption holds then the following result can be proved, by means 
of a slight adaptation of the results presented in j8j. 

Proposition 3 Let (A k ,B k ,C k ) e R m * xm fc x R m fe xl x R lxmfc be a triplet 
of matrices in prime form. Furthermore, let g k '■ R mfc — > R a bounded 
locally Lipschitz function that agrees with g k (-) on Br = {£ € R mfe : |£| < 
R}, let Di = diag(^,^ 2 , . . . ,£ mk ) with £ a positive design parameter, and let 
(c , . . . , c mk -i) be such that the polynomial A mfe + c A mfc_1 + . . . + 
Hurwitz. Then there exist R > and £* > such that for all £ > £* the 
triplet (ip k (-),j k (-),G k ) defined as 

ipk(i]k) = A k rj k + B k g k (n k ) , 7fc(??fc) = C k r] k , 

G k = Di coI(cq, . . . , c mk -i) has the asymptotic internal model property. 



5 Simulation results 

In this section we present simulation results about the theory presented in 
Section III by considering two different nonlinear oscillators as system ([6]), 
namely a Van der Pol and a Duffing oscillator. We consider the case of three 
agents (N = 3). 

In case of Van der Pol, system ©, $7$) takes the form 



The topology of the simulated network is described by the incidence matrix 



with the eigenvalues of the corresponding Laplacian that fulfill (|12p with 
[i = 1.4. Figures [TJ [2] present the simulation results obtained with the three 
systems of the form (|14p where the coupling term has been fixed as in Lemma 
1 with g = 10 and a = 1. In particular, Figure Q] shows the phase-plane of 
the three oscillators initialized respectively at wi = (1,1), W2 = (3,3) and 
W3 = (5,5), while Figure [2] shows the time behaviors of their first state 
variable. 



Wi = W2 

W2 = 2(1 — w\)w2 — W\ 



y* = w\. 





Figure 1: Phase-plane of the three networked Van der Pol Oscillators. 



Figure 2: Time behavior of the first state variable of the three Van der Pol 
Oscillators. 



As a second example we consider the case of a Duffing oscillator described 

by 

w 2 = -2wi +wf y 1 ' 

The topology of the network is described by the same incidence matrix 
considered in the previous example with the same parameters that have 
been used in the coupling terms. Figure [3] shows the phase-plane of the three 
oscillators initialized respectively at w\ = (1, 1), w 2 = (3, 3) and W3 = (5, 5), 
while Figure U] shows the time behaviors of their first state variable. 



6 Conclusions 



The problem of reaching a consensus among the outputs of a set of networked 
nonlinear systems was considered. The output reference signal is thought 
of as generated by the autonomous nonlinear exosystem of the form ([5]). A 
first result of the paper is presented in Section III where it is shown how 
TV diffusively-coupled exosystems of the form (I10p reach a consensus over 
a trajectory of ([5|). This result is instrumental to the theory presented in 
Section IV where the nonlinear output regulation theory is adopted to design 
local regulators that make the outputs of N heterogenous systems to track 
a common reference signals solution of ([5]). 



Figure 3: Phase-plane of the three networked Duffing Oscillators. 



Figure 4: Time behavior of the first state variable of the three Duffing 
Oscillators. 
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